On the branching ratio of the "second class" r — > rj'TTiyr decay 
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Abstract 

We present an assessment of the standard model expectations for the branching ratio of the 
isotopic spin and G-parity violating decay r — )• ri'irvr- This estimate is based on a vector and 
scalar meson dominance parametrization of the relevant form factors, that explicitly accounts for 
TT^ — rj — rj' mixing. The numerical results obtained in this framework indicate a branching ratio 
one order of magnitude (or more) below the current experimental limit, and suggest the possibility 
of evidencing some novel interaction in high statistics studies of this decay. 
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The semileptonic transition r — r^'Trz/,- belongs to the category of the so-called "second 
class current" decays This kind of processes can in the standard model occur only 
via an isotopic spin and G-parity violation mechanism and, consequently, are strongly sup- 
pressed by the small d — u current quark mass difference. The present experimental limit, 
BR(r — )■ ti't^Vt) < 7.2 x 10~^ {2!, is rather constraining. Nevertheless, depending on the ac- 
tual smallness of the standard model branching ratio, this upper limit might still be high 
enough to allow some room available for revealing "nonstandard" effects in high statistics 
studies of this decay, for example, at a Super B factory j^. This has recently been empha- 
sized with regard to r — )■ r/Trz/^ in Refs. , and earlier in Refs. |6|, l7|, and should be the 
case also of the decay of interest here. 

While there are several theoretical estimates of r — )■ tjttUt-, based either on the direct 



application of p(770) and ao(980) pole dominance 



J, |8| or on current algebra and chiral 



perturbation theory j^, [lo|, recent attempts to estimate the branching ratio for r — > ri'iru^ 
appeared in the literature are, to our knowledge, much fewer [Uj. In this note, we extend 
to this process the method used for the calculation of r — )■ rj-KUr in Ref. This approach 
consists of a parametrization of the spin-1 and spin-0 relevant form factors in the timelike 
region, in terms of p(770),p' = p(1450) and ao(980),aQ = ao(1450) exchanges, respectively, 
with coupling constants either theoretically estimated within the quark model or limited 
from the available phenomenology, and constrained by the value of the isotopic spin violating 
— 7] mixing parameter earlier evaluated from chiral symmetry breaking. 
Indeed, in this regard, the scale of the transition amplitude for r — t- rj-KUr is set by a 7r° — 77 
mixing parameter calculable in chiral perturbation theory to leading and next-to-leading 



order, see Refs. 



12 



13| and the application in Ref. {q]. The scale of r — )■ r^'Trz/,- can be set by 



a 7r° — ?7 — 77' mixing parameter e,,/^ that, the t]' not being a Nambu-Goldstone boson except 
for Nc ^ 00, can phenomenologically be estimated in a chiral symmetry breaking scheme 
supplemented by the determination of the t] — rj' mixing angle [l^. The other, obvious, 
difference between the two processes is a kinematical one, namely, in the case of r — )■ ri'irur 
the 1~ and 0^ exchanged ground states are either well-below or just around (if the width is 
taken into account) the threshold of the decay phase space. 

With = u'y^d the weak vector current, the hadronic matrix element for r — )■ ri'irUj. can 
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be decomposed into spin-1 and spin-0 exchange form factors, respectively, as 



{7r^ikWip)\V,\0) = -V2 



— + fo{t)- 



(1) 



t ■ ""^ ' t 

Here, q = p + k and t = is the invariant mass squared of the emitted rj'ir pair. 

With to = {Mrj' + Mtt)^ and A(x, 2;) = + + — 2{xy + yz + zx), the partial width 

is: 



r(r+ ^ r]'n+Ur) 



t3 



X 0/iWlM(2t + M,2)A(t,M^2,,M,2)) + |/o(t)p3M,^(M^2,-M,2)2] (2) 
For the spin-1 form factor we assume, along the lines of Refs. fisl . the unsubtracted 



linear combination of p and p' polar forms 



flit) 



f p9pri'-K 



M. 



^'Ml,-t-iMp,Tp.{t) 



(3) 



17| 



(4) 



Here, denoting by M and F the mass and total width of a spin-L resonance 

where in our case q is the momentum in the ri'ir cm. frame and L = 1,0 for spin-1 and 
0, respectively. Notice that, in Eq. ([3]), we have omitted the p width, this pole is below 
the threshold to and the numerical results are almost insensitive to this approximation. 
In Eq. (|3]), fp and (yfp^/^ are the vector meson couplings to e"'"e~ and to 77'??, respectively, 
and the coefficient [3p = {Mp/MpiY x {fp'gp'rj'Tr/ fpOpri'Tr) parametrizes the contribution of the 
radial excitation p'. In some sense, Eq. ([3]) resembles the modification of the p propagator 
introduced in Ref. [isl ]. 

We incorporate isotopic spin violation through the rj' -H- 7r° transition (7r''|H'|?7'). In 
the quark model, neglecting isospin breaking of electromagnetic origin, H' = -AmS^ with 
Am = rrid — rriu and the scalar densities 5** = q^q (A's are Gell-Mann matrices). To leading 
order in Am, the {pri'ir) and (p'tj'tt) trilinear couplings will be assumed to be proportional 
to the {pTcir) and (p'tttc) ones, through the mixing parameter 

(7rO|H'|r7') 



(5) 



so that 



dprj'-K ^q' -kQ pn-K ) 



Qp'ri'-K ^q' -kQ p' TTTT- 



(6) 



In the soft-pion limit, by applying current algebra: 

Am 



{0\V2Po + P8W), 



(7) 



with Fj, the pion decay constant and the pseudoscalar densities P* = q\^^q. Similar 
relations as in Eqs. ([5]) and ([7]) can be derived also for the case of the rj. 

For the rj and rj' states we adopt the single mixing angle scheme in terms of the pure 
octet and singlet pseudoscalars, 



\rj) = cos^It^s) — sin 9\rio); \r]') = sin^lr^g) + cos 6*1770). 



(8) 



By relating the P* in Eq. ([7]) to the divergences of the octet and singlet axial currents (f/(l) 
anomaly included), in the SU{2) x SU{2) limit, after some algebra one obtains the relation 



{0\V2Po+PsW) 



^Pg^ [cos^ OM^ + sin^ OM^, - 2 cos^ 6 (M^, - M^)] . (9) 



We take for the mixing angle the value 6 = —20° 19|]. Values of 6 in the range, say, 
[—15°, —20°] have been obtained recently 20|, depending on the physical process and on the 
model, and also two-angles rj — rj' mixing schemes have been proposed [21j]. On the other 
hand, the dependence of Eq. ([9]) on 6* is rather mild. By combining Eqs. (I5])-(I9]), assuming 
the quark mass ratios rriu/md ~ 0.55 and ms/md ~ 18.9 [22|, and the pseudoscalar decay 
constants Pg ~ we find the approximate value en'-n — 3 x 10^'^. To encompass the other 



determinations 
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23 



and somehow account for theoretical uncertainties, we in the 



sequel will allow for e^/^r the range of values 

(3± 1) X 10^1 



(10) 



Notice that, with these values, and (^ptttt = 6 from the experimental p- width, the determi- 
nation of the [prj'Tc) coupling constant from Eq. (jG]) is somewhat smaller than obtained in 



Ref. 



25|, and falls well-below the upper limit derived in Ref. 11|. 



As regards the estimate of the parameter /3p in Eq. ([3]), a calculation within the con- 
stituent quark model, using the ratio of p and p' wave functions at the origin 26|, indicates 
fpi/fp ^ 1.1. Moreover, an assessment (actually, an upper limit) of gpi-^T, can be obtained 
by identifying V{p' — t- tttt) to V pi ~ 0.4 GeV. This gives for [3p the upper limit: |/3p| < 0.18. 
We finally complete the numerical input needed for the parametrization ([3]) by taking, from 



experimental data on p decays, fpgpT,^/Mp ~ 1.2 [2 



We turn to the spin-0 exchange form factor fo{t), determined by the matrix element of 
the divergence of the vector current in ([1]), through the relation exhibiting the isotopic spin 
violation proportional to Am: 



We assume a similar parametrization as in Eq. (|3]), dominated by oq s'^nd Og poles: 



Mt) 



FaoQaori'iT 

Ml, - Ml 



M2 



Ml-t-lMa,Ta, 



+ /5a 



M2 



M"^, - t- iMa' Ta' U) 



(12) 



Here, the factor f3a multiplying the Qq pole is defined by the ratio: (3a = Fa'^ga'^ri'n/ Faogaor]'-K, 
with gaor]'n and (^a^^'Tr the trilinear (aor/'Tr) and (aQr^'vr) coupling constants, respectively, and 

- V2Fa,Ml = {ao{q)\zd^V,\0), -v^F.^M^ = {a',{q)\td'^V,\0) . (13) 

In Eq. ( IT2I) . we neglect the t- variation of in order not to generate, with Ma^ so close to 
the threshold to 5 spurious imaginary parts. 

In a U{3) scheme with rj — rj' mixing according to ([8]), one obtains for 6 

(V2cos6' + sin 6* 



-20°: 



9aorin 



daLriTT 



(cos 9 — \/2si 



0.70. 



sm I 



(14) 



This relation is in agreement with the experimental ratio T{aQ — > ?7'7r)/r(ao — ?■ r^vr) ~ 

Along the lines of Ref. fs], to exploit Eq. f[T^ we adopt the 



0.35 121, see also Ref. 



experimental value gaorjiT — 2.80 GeV accompanied by a total width Tag — 100 MeV |29|. 
For the ga'^rin coupling constant, we assume that the t]tt, tj'tt, KK and wvrvr decay channels 
of the Oq saturate the total width Ta' = 265 MeV. Using the ratios between partial widths 



T{ri'TT)/T{riTT), T{KK)/T{r]TT) and r(a;7r7r)/r(?77r) reported in [27|, one would find r(aQ — >• 
rjn) ~ 20 MeV, and the consequent estimate ga'rjw — 1-32 GeV. Notice that these values 



agree to a good extent with the predictions recently derived in Ref. (30| from an SU(6) 
breaking approach, and in Ref. |3l| using QCD sum rules. 

To estimate f3a, we further need the values of the the constants and Fa' in ([T3|) . 
and we assume both 0^ scalars to be p-wave qq states. Similar to [11], one might adopt 
the SU (6) framework and relate, by means of current algebra equal-time commutators and 
single-particle saturation of the ensuing sum rules, Fa^ to the analogous l"*" constant Fa^. 



In the SU{2) x SU{2) limit, the 2nd Weinberg sum rule [32| would then relate Fa^ = Fp. 



Using the same procedure for F^^, we would obtain the ratio 

Fa' Ml F, 



_^ ^0 



ations 



with Fpi /Fp previously estimated, hence the estimate = 0.23. 

Anther possibility would be to use the QCD sum rules and local hadron duality re 
M|(n) = {n + 1)M| and = F|/(n + 1), where S denotes 0+ scalar |33|, b^]. For 

n = 1, the former relation is in perfect agreement with the experimental and iVf^j,; the 
latter one indicates Fa'^/Fa^ = 0.70. Numerically, with this approach we get the estimate 
IPal — 0.33, which we ultimately adopt for the final numerical calculations and, anyway, is 
compatible within our uncertainties to the value obtained before. 

Finally, theoretical estimates of Fag are available, and we take the values Fa^ ~ 1.3 — 



1.6 MeV obtained from QCD sum rules 35 1. 



Inserting in Eq. (|2]) the parametrizations (|3]) and f lT2|) with the input coupling constants 



derived above, in particular with e^/^ as in (ITU]) . \/3p\ < 0.18, consistent with [36|, and 
\Pa\ < 0.33, we obtain predictions for the decay rates from spin-1 and spin-0 exchanges. For 
any fixed pair of f3p, f3a in the mentioned ranges, we fix the corresponding value of e^/7r lying 
in the range ( !T0|) . such that, as needed to cancel the kinematical singularity in Eq. ([1]), the 
condition /i(0) = /o(0) is numerically satisfied. From this procedure, we find the following 
intervals for the spin-1 and spin-0 exchanges: 

BR(r ^ r/'7ri/^)|,pin-i = 1.4 x 10"^ - 3.4 x 10"^ (16) 

and 

BR(r ^ ri'ni^r)\spm~o = 6.0 x 10"® - 1.8 x 10"^. (17) 

With the current, limited, knowledge of the p' and Oq properties, and the generous limits 
allowed to (3p and we are not able to predict the branching ratio of r — )• r(Txv^, with a 
fully theoretical calculation, to a better accuracy. As one can see, Eqs. (fT6l) and (fT7|) clearly 
reproduce the general expectation of the phase space suppression of the spin-1 exchange vs. 
the spin-0 one. The upper values in these equations fall below the upper limit of the order 
of 10~^ presented in ll|, and should encourage high statistics studies of this process in the 
quest for "non-standard" exchanges. 
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